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Abstract. For a discrete mechanical system on a Lie group G determined 
by a (reduced) Lagrangian £ we define a Poisson structure via the pull-back 
of the Lie-Poisson structure on the dual of the Lie algebra g* by the corre- 
sponding Legendre transform. The main result shown in this paper is that 
this structure coincides with the reduction under the symmetry group G of 
the canonical discrete Lagrange 2-form oil on G x G. Its symplectic leaves 
then become dynamically invariant manifolds for the reduced discrete system. 
Links between our approach and that of groupoids and algcbroids as well as 
the reduced Hamilton-Jacobi equation are made. The rigid body is discussed 
as an example. 
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1. Introduction 

Background. This paper continues our development of discrete Lagrangian me- 



chanics on a Lie group introduced in Marsden, Pekarsky, and ShkoUer |MPeS 



In our earlier paper, using the context of the Veselov method for discrete mechan- 
ics, discrete analogues of Euler-Poincare and Lie-Poisson reduction theory (see, for 
example, Marsden and Ratiu [ MR 99 ) were developed for systems on finite dimen- 
sional Lie groups G with Lagrangians L : TG E that are G-invariant. The 
resulting discrete equations provide "reduced" numerical algorithms which mani- 
festly preserve the symplectic structure. The manifold G x G is used as the discrete 
approximation of TG, and a discrete Lagrangian L : G x G ^ R is constructed from 
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a given Lagrangian L in such a way that the G-invariance property is preserved. 
Reduction by G results in a new "variational" principle for the reduced Lagrangian 
£ : G = (G X G)/G — >• R, which then determines the discrete Euler-Poincare (DEP) 
equations. Reconstruction of these equations is consistent with the usual Veselov 
discrete Euler-Lagrange equations developed in |WM 97 , MPS 98 1, which are nat- 



urally symplectic-momentum algorithms. Furthermore, the solution of the DEP 
algorithm leads directly to a discrete Lie-Poisson (DLP) algorithm. For example, 
when G — S0(7i), the DEP and DLP algorithms for a particular choice of the dis- 



crete Lagrangian L are equivalent to the Moser- Veselov [MoV 91] scheme for the 
generalized rigid body. 

Main Results of this Paper. We show that when a discrete Lagrangian L : G x 
G ^ M is G-invariant, a Poisson structure on (a subset) of one copy of the Lie group 
G can be defined which governs the corresponding discrete reduced dynamics. The 
symplectic leaves of this structure become dynamically invariant manifolds which 
are manifestly preserved under the structure preserving discrete Eulcr-Poincare 
algorithm (see Section 2.1). 

Moreover, starting with a discrete Euler-Poincare system on G one can readily 
recover, by means of the Legendre transformation, the corresponding Lie-Poisson 
Hamilton- Jacobi system on g* analyzed by Ge and Marsden |]GM 88|| ; the relation- 
ship between the discrete Euler-Lagrange and discrete Euler-Poincare equations and 
the Lie-Poisson Hamilton- Jacobi equations was examined from a different point of 



view in our companion paper |MPeS 99 1 



We also apply Weinstein's results on Lagrangian mechanics on groupoids and 



algebroids | W 96 to the setting of regular Lie groups. The groupoid-algebroid 
setting reveals new and interesting connections between discrete and continuous 
dynamics. 

2. Discrete reduction 
In this section we review the discrete Euler-Poi ncare reduction of a Lagrangian 



system on G x G considered in detail in | MPeS 99 1. We approximate TG by G x G 



and form a discrete Lagrangian L : G x G ^ M from the original Lagrangian 
L : TG ^ R by 

^{9k,9k+i) = L{K{gk, gk+i),x{9k, gk+i)), 

where k and x ^re functions of {gk, gk+i) which approximate the current configura- 
tion g{t) G G and the corresponding velocity g{t) G TgG. We choose discretization 
schemes for which the discrete Lagrangian L inherits the symmetries of the original 
Lagrangian L: L is G-invariant on G x G whenever L is G-invariant on TG. In 
particular, the induced right (left) lifted action of G onto TG corresponds to the 
diagonal right (left) action of G on G x G. 

Having specified the discrete Lagrangian, we form the action sum 

N~l 

§ = X! H9k,gk+i), 

k=0 

which approximates the action integral S = J Ldt, and obtain the discrete Euler- 
Lagrange (DEL) equations 

D2U9k^i,9k) + DM9k,gk+i) = (2.1) 
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as well as the discrete symplectic form wl, given in coordinates on G x G by 



dgldgi 



k+l 



In (2.1), Di and D2 denote derivatives with respect to the first and second argu- 
ment, respectively. The algorithm (2.1) as well as wl are obtained by extremizing 
the action sum S : G^'^^ K with arbitrary variations. Using this variational 
point of view, it is known that the flow Ft of the DEL equations preserves this 
discrete symplectic structure. This result was obtained using a discrete Legendre 
transform and a direct computation in | V 88|, V 91, WM 97] and a proof using the 



variational structure directly was given in |MPS 98 



Remark 2.1. We remark that the discrete symplectic structure is not globally 
defined, but rather need only be nondegenerate in a neighborhood of the diagonal 
A in G X G, i.e., whenever g^ and gk+i are nearby. Section 3 of |MPS 98 1 shows that 



wl arises from the boundary terms of the discrete action sum restricted to the space 
of solutions of the discrete Euler-Lagrange equations; an implicit function theorem 
argument relying on the regularity of the discrete Lagrangian L is required in order 
to obtain solutions to the discrete Euler-Lagrange equations, and this regularity 
need only hold in a neighborhood of the diagonal A C G x G. 

2.1. The discrete Euler-Poincare algorithm. The discrete reduction of a right- 
invariant system proceeds as follows (see [MPeS 99 for details). The case of left 
invariant systems is similar. Of course, some systems such as the rigid body are 
left invariant. 

The induced group action on G x G by an element ^ G G is simply right multi- 
plication in each component: 

9 ■ {9k,gk+i) 1-^ {gkg,gk+ig), 

for all gk,gk+i e G. 

The quotient map is given by 



7rd:GxG^(GxG)/G-G, (5fc, ^ (2.2) 

One may alternatively use gk+ig^^ instead of gkg^+i as the quotient map; the 
projection map (2.2) defines the reduced discrete Lagrangian £ : G ^ M for any 
G-invariant L by ^ o tt^j = L, so that 

Hdkgkli) = Hgk,gk+i), 

and the reduced action sum is given by 

s = J2 Kfkk+i), 

k=0 

where fkk+i = gkg^+i denotes points in the quotient space. A reduction of the 
DEL equations results in the discrete Euler-Poincare (DEP) equations 

Rl.,/(fkk+i)-L)^_J{fk-ik)^0 (2.3) 

for k — l,...,Af — 1, where R*^ and L*^ for f E G are the right and left pull- 
backs by /, respectively, defined as follows: for ag G T*G, R*fag G q* is given 
by {R*jag,£) = {ag,TRf ■ ^) for any ^ € g, where TRf is the tangent map of the 
right translation map Rf : G G; h ~* hf, with a similar definition for L^-. Also, 
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G —> T*G is the differential of I defined as follows. Let g'^ be a smooth curve in 
G such that g" — g and {d/de)\e=og'^ = v. Then 

e'{g)-v={d/de)U^oi{9')- 



For the other choice of the quotient in (2.2) given by /ifc+ifc = gu+idk ^ ^ the DEP 
equations are 

Lk+J'i^k+ik) - Rl^^_/{hkk-i) = (2.4) 

Remark 2.2. In the case that L is left invariant, the discrete Euler-Poincare equa- 
tions take the form 

^}.+..^'(^+i^)-^L-/'(^^"i)=0 (2.5) 

where fk+ik = 9k+i9k is in the left quotient (G x G)/G. 

Notice that equations (2.4) and (2.5) are formally the same. 

We may associate to any G"^ function F defined on a neighborhood V of A C 
G X G its Hamiltonian vector field Xp on V D A satisfying Xp-^oji^ = dF , where c?F, 
the differential of F, is a one-form. The symplectic structure ljl naturally defines 
a Poisson structure on a neighborhood V of A (which we shall denote {•, -Igxg) by 
the usual relation 

{F, H}gxg = ^^l{Xf, Xh). 
Theorem 2.2 of MPeS 99[ | states that if the action of G on G x G is proper, the 
algorithm on G defined by the discrete Euler-Poincare equations (2.3) preserves the 
induced Poisson structure {•, -jc on W C G given by 

{f-.hjcond^ {f oTTd,hond}GxG (2.6) 

for any G^ fmictions /, h onU, where U = TTd{V). 

Using the definition fkk+i ~ 9kgk+i; the DEL algorithm can be reconstructed 
from the DEP algorithm by 

(ff/c-i, 5fc) ^ {9k,gk+i) = {.fk-ik ■ 9k-iJkk+i ■ gk), (2.7) 

where fkk+i is the solution of ( ^.3[ ) . Indeed, f^7|}_^_l -gk is precisely gk+i- Similarly one 
shows that in the case of a left G action, the reconstruction of the DEP equations 
( p^ ) is given by 

igk-i,gk) {9k,gk+i) = (gk-i ■ fkk-i,9k ■ fk+ik)- 

2.2. The discrete Lie-Poisson algorithm. In addition to reconstructing the 
dynamics on G x G, one may use the coadjoint action to form a discrete Lie- 
Poisson algorithm approximating the dynamics on q* | MPeS 99| 

Mfc+i-Ad}^^^^.Mfe, (2.8) 

where ^fc := Ad*-i fiQ is an element of the dual of the Lie algebra, /zq is the constant 

of motion (the momentum map value), and the sequence {fkk+i\ is provided by 
the DEP algorithm on G. 

The corresponding discrete Lie-Poisson equation for the left invariant system is 
given by 

Hfc+i-Ad* 1 -Hfe, (2.9) 

J k + lk 

where 11^ := Ad*^ ttq and ttq is the constant momentum map value. Henceforth, we 
shall use the notation € 0* for the right invariant system and 11 S fl* for the left. 
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3. POISSON STRUCTURE AND INVARIANT MANIFOLDS ON LlE GROUPS 

Discretization of an Eulcr-Poincare system on TG results in a system on G x G 
defined by a Lagrangian L. If it is regular, the Legendre transformation (in the 
sense of Veselov) Fh define a symplectic form (and, hence, a Poisson structure) 
on V C G X G via the pull-back of the canonical form from T*G. Then, general 
Poisson reduction applied to these discrete settings defines a Poisson structure on 
the reduced space U = Trd{V) C G. This approach was adopted in Theorem 2.2 of 
|MPeS 9^ . 



Alternatively, without appealing to the reduction procedure, a Poisson struc- 
ture on a Lie group can be defined using ideas of Weinstein | W 96| on Lagrangian 



mechanics on groupoids and their algebroids. The key idea can be summarized in 
the following statements. A smooth function on a groupoid defines a natural (Le- 
gendre type) transformation between the groupoid and the dual of its algebroid. 
This transformation can be used to pull back a canonical Poisson structure from 
the dual of the algebroid, provided the regularity conditions are satisfied. 

The ideas outlined in this section can be easily expressed using the groupoid- 
algebroid formalism. Such a formalism is suited to the discrete gauge field theory 
generalization as well as to discrete semidirect product theory; nevertheless, the 
theory of groupoids and algebroids is not essential for the derivations, but rather 
contributes nicely to the elegance of the exposition. 

3.1. Dynamics on groupoids and algebroids. In this subsection, we show 
that our discrete reduction methodology is consistent with Weinstein's groupoid- 
algebroid construction; the contents of this subsection are not essential for the 
remainder of the paper. 



We briefly summarize results from Weinstein |W 96| and refer the reader to the 
original paper for details of proofs and definitions. Let F be a groupoid over a set 
M, with a, /3 : r — > M being its source and target maps, with a multiplication map 
m : r2 ^ r, where r2 = {{g, /i) G F x F | /3(g) = a{h)}. Denote its corresponding 
algebroid by A. 

The Lie groupoids relevant to our exposition arc the Cartesian product G x G 
of a Lie group G, with multiplication {g,h){h,k) = {g,k), and the group G itself. 
The corresponding algebroids are the tangent bundle TG and the Lie algebra g, 
respectively. The dual bundle to a Lie algebroid carries a natural Poisson structure. 
This is the Poisson bracket associated to the canonical symplectic form on T*G and 
the Lie-Poisson structure on g*, respectively. 

Lagrangian mechanics on a groupoid F is defined as follows. Let £ be a smooth, 
real-valued function on F. £2 the restriction to F2 of the function (g. h) t-^ C{g) + 
C{h). 

Definition 3.1. Let E£ C F2 be the set of critical points of C2 along the fibers of 
the multiplication map m; i.e. the points in Y,c stationary points of the function 
C{g) + C(h) when g and h are restricted to admissible pairs with the constraint that 



the product gh is fixed |W 96| 



A solution of the Lagrange equations for the Lagrangian function £ is a 
sequence . . . , 5-2, 50j 5ii 52, ■ ■ • of elements of F, defined on some "interval" in 
Z, such that {gj,gj^i) S E£ for each j. 

The Hamiltonian formalism for discrete Lagrangian systems is based on the fact 
that each Lagrangian submanifold of a symplectic groupoid determines a Poisson 
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automorphism on the base Poisson manifold. RecaU that the cotangent bundle T*r 
is, in addition to being a symplectic manifold, a groupoid itself, the base being A* ; 
notice that both manifolds are naturally Poisson. The source and target mappings 
a, P : T*T A* are induced by a and 

Definition 3.2. Given any smooth function C on T, a Poisson map Ac from A* 
to itself which may be said to be generated by C is defined by the Lagrangian 



submanifold dC(T) (under a suitable hypothesis of nondegeneracy) [W 96| 



The appropriate "Legendre transformation" FC in the groupoid context is given 
hy a o dC : r ^ A* oi f3 o dC : T ^ A* , depending on whether we consider right or 
left invariance (through the definition of maps a and The transformation FC 
relates the mapping on F defined by E£ with the mapping on A* . FC also pulls 
back the Poisson structure from A* to F, which, in general, is defined only locally 
on some neighborhood U C T. In the context of a Lie group, this means that any 
regular function £ : G ^ R defines a Poisson structure on U. We shall address this 



issue in the next subsections. The reader is referred to |W 96 1 for an application of 
the above ideas to the groupoid M x M when the manifold M does not necessarily 
have group structure. 

3.2. DEP equations as generators of Lie-Poisson Hamilton- Jacobi equa- 
tions. A Lie group G is the simplest example of a groupoid with the base being 
just a point. Its algebroid is the corresponding Lie algebra g, with the dual being 
Q* . Consider left invariance and let a general function C on the group be specified 
by the discrete reduced Lagrangian (. : G ^ M.. Then, the Legendre transform 
defined above is given by 

F£ = L*gOd£ -.G^s*, 
where d£ : G ^ T*G. Using these transformations we define 

nfe_i = F£{fkk-i) = L%^_^od£ifkk-i)- 
Recall the DEP equation ( p.5D for left-invariant systems : 

L%^.J^ih+ik)-R%,_,d£ifkk^,)^0, 

where we have identified the notations £' and d£. The latter equation can be 
rewritten as a system 

nk = L}od£{f), 

Uk+i=R}od£if), ^''■'> 

where the first equation is to be solved for / (which stands for fk+ik) which then 
is substituted into the second equation to compute Ilfc+i. 

This system is precisely the Lie-Poisson Hamilton- Jacobi system described in 



GM 88[ with the reduced discrete Lagrangian £ playing the role of the generating 
function. This means that there is no need to find an approximate solution of the 
reduced Hamilton- Jacobi equation |GM 88 1. Notice also that the DLP equation 



(2.£) is a direct consequence of the system (3.1) 



Hfc+i==Ad* 1 -Hfe. 

J k+lk 
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The following diagrams relate the dynamics on G and on g*: 



G 



G 



Ft 



fkk-l 
Fl 



0' 



0' 



n, 



fk+lk 
Ft 



(3.2) 



where and are given in Definitions 3.1 and 3.2 



3.3. Some Advantages of Structure-preserving Integrators. As we men- 
tioned above, the "Legendre transform" F£ allows us to put a Poisson structure 
on the Lie group G, which, of course, depends on the discrete Lagrangian L on 
G xG, and hence on the original Lagrangian L on TG (if we consider this from the 
discrete reduction point of view) . It follows that the reduction of the discrete Euler- 
Lagrange dynamics on G x G is necessarily restricted to the symplectic leaves of 
this Poisson structure, so that these leaves are invariant manifolds, and correspond 
(under Fi) to the symplectic leaves (coadjoint orbits) of the continuous reduced 
system on g*. 

These ideas are the content of the following theorems. 

Theorem 3.1. Let L be a right invariant Lagrangian on TG and let L be the 
Lagrangian of the corresponding discrete system on V C G x G. Assume that L is 
regular, in the sense that the Legendre transformation _FL : V — > -FL(V) C T*G is 
a local diffeomorphism, and let the quotient maps he given by 

TTd-G xG ^ {G xG)/G^G and ir : T*G ^ {T*G)/G q* . 

Let i be the reduced Lagrangian on G defined by 

L = £ O TTd, 

and let 

Fe : U C G ^ g* 

be the corresponding Legendre transform. Then the following diagram commutes: 

FL 



VcGxG 
U cG 



T*G 



(3.3) 



Fe 



Proof. First, we choose coordinate systems on each space. Let {gk,gk+i) G G x 
G and {g,p) G T*G, so that the discrete quotient map ( |2.2| ) is given by tt^ : 
{gk,gk+i) ^ fkk+i = 9k9klii and the continuous quotient map by tt : {g,p) ^ 
fi = RgP. Recal l that th e fiber derivative FL in these coordinates has the following 
form (see, e.g., [IWM 97| ) 



FL:GxG-^r*G; {gk,gk+i 
Then the above diagram is given by: 



{gk,Di'L{gk,gk+i))- 



(3fc,.g/c+i) 



{gk,Pk 



dgk' 



(3.4) 



^„-i 9k 



M = R*g,Pk, 
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where / stands for fkk+i — QkQ^^i- To close this diagram and to verify the arrow 
determined by F^ compute the derivative of L using the chain rule 

where we have used that, according to the definition of /, the partial derivative 



— — is given by the linear operator TR -i . (3.5) is precisely the Legendre trans- 
ogk ^'==1 

formation Ft for a right invariant system (see the previous subsection). □ 

Corollary 3.1. Reconstruction of the discrete Lie-Poisson (DLP) dynamics on g* 
by TT^^ corresponds to the image of the discrete Euler- Lagrange (DEL) dynamics 
on G X G under the Legendre transformations Fh and results in an algorithm on 
T*G approximating the continuous flow of the corresponding Hamiltonian system. 

Proof. The proof follows from the results of the previous subsection, in particular. 



diagram (3.2) relates the DLP dynamics on g* with the DEP dynamics onU C G 



which, in turn, is related to the DEL dynamics on V C G x G via the reconstruction 



(2.7). □ 



An important remark to this corollary which follows from the results in [ KMO 99| 



(see also |KMOW 99) |) is that, in general, to get a corresponding algorithm on the 
Hamiltonian side which is consistent with the corresponding continuous Hamilton- 
ian system on T*G, one must use the time step /i-dependent Legendre transform 
given by the map 

(ff/c,5fc+i) ^ {gk,-hDih{gk,gk+i))- 

The results of this paper are not effected, however, as we assume h to be constant 
and so we would simply add a constant multiplier to the corresponding symplectic 
and Poisson structures. For variable time-stepping algorithms, this remark is crucial 
and must be taken into account. 

Theorem 3.2. The Poisson structure on the Lie group G obtained by reduction of 
the Lagrange symplectic form wl on V C G X G via tt^ coincides with the Poisson 
structure onlA <Z G obtained by the pull-hack of the Lie-Poisson structure uj^ on q* 
by the Legendre transformation F£ (see diagram ( fg. j[ i above). 



Proof. The proof is based on the commutativity of diagram (|3.3| ) and the G invari- 
ance of the unreduced symplectic forms. Notice that G and g* in ( |3.3| ) are Poisson 
manifolds, each being foliated by symplectic leaves, which we denote E / and for 
/ € G and /i g g*, respectively. Denote by cj/ and the corresponding symplectic 
forms on these leave s. B elow we shall prove the compatibility of these structures 
under the diagram ( |3.3| ). Repeating this proof leaf by leaf establishes then the 
equivalence of the Poisson structures and proves the theorem. 

Recall that the Lagrange 2-form on V C G x G derived from the variational 
principle coincides with the pull-back of the canonical 2-form LUca.n on T*G (see. 



e.g., I MPS 98, WM 97|). Recall also that for a right-invariant system, reduction of 
T*G to g* is given by right translation to the identity e G G, i.e. any p G T*G is 
mapped to fi — R*p G g* = T*G. Thus, for any g e 7r~^(/i), where ^ e g*. 
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SO that (tt ^)* — {R*_i)* puUs back Wcan to w^. Henceforth, tt ^ shall denote the 
inverse map of tt restricted to TgG* . 

Let us write down using the above notations how the symplectic forms are being 
mapped under the transformations in diagram (3^); we see that 



V c G X G 



T*G 



UcG 



Fl 



Ft 



^can 



(3.6) 



Then, using the coordinate notations of diagram (3.4), for any f £ G and u,w G 
U}{f){u,v) = LJ^,{^i){TFl{u),TF^(v)), (3.7) 



where /i = Ft{f) G g*. Continuing this equation using diagram (B.C), we have that 
w/(/)(u, v) - c^can((.gfc,Pfc))(T7r-i o TFi{u),TTT-^ o rF^(t;)) 

= c^L((5fe, 5fe+i))(rFL-i „ y^-i ^ y^^(y)^ j^^L-i o Ttt-i o T^^£(t;)), (3.8) 

where (gk,Pk) G ""^^(z^) a-nd Ttt^^ denotes TR*_i. 
Using (^), it follows that 

F£oTTd = TroFh 
and, hence, for the tangent maps, we have that 

TFe o TTTrf = Ttt o Ti^L. 



So, if u, u in (3/7) are images of some G i nvar iant vector fields U,V onV C G x G, 
i.e. u = TTTd{U),v = Tnd{V), then from (|U) it follows that 

^/(/)(w,w) = tJL((5'fc,.gfc+i))(t/, V"), 

where {gk,gk+i) = '^d^if) ^^'^ U,V £ Ti^g^g^^^-^G x G. The last equation precisely 
means that ojf is the discretely reduced symplectic form, i.e. the image of cjl under 
the quotient map -Kd- □ 

Analogous theorems hold for the case of left invariant systems. 

More General Configuration Spaces. Similar ideas carry over to the integra- 
tion of systems defined on a general configuration space M with some symmetry 
group G. In this case, the reduced discrete space (M x M)/G inherits a Poisson 
structure from the one defined on M x M (analogously to ( |2.6| )). Its symplectic 
leaves again become dynamically invariant manifolds for structure-preserving inte- 
grators and can be viewed as images of the symplectic leaves of the reduced Poisson 
manifold T*M/G under appropriately defined "Legendre transformations". This is 
a topic of ongoing research that builds on recent progress in Lagrangian reduction 
theory; see [ |MRS 99t . 

3.4. Poisson structures of the rigid body. As an example of applications of 
the above ideas, we consider the dynamics of the rigid body and its associated 
reduction and discretization (see, e.g. jMR 99| , [MoV 91| , |LS 96| , [MPeS 99[ for more 
details). 

The Basic Set Up. The configuration space of the system is S0(3). The cor- 
responding Lagrangian is determined by a symmetric positive definite operator 
J : 5o(3) so(3), defined by J(^) = A^ -I- ^A, where ^ G so(3) and A is a diagonal 
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matrix satisfying A.^ + Aj > for all i ^ j. The left invariant metric on SO (3) is 
obtained by left translating the bilinear form at e G SO (3) given by 

(c,o = jTr(e^j(e)). 

The operator J, viewed as a mapping J : so(3) — > so(3)*, has the usual interpre- 
tation of the inertia tensor, and the A^ correspond to the sums of certain principal 
moments of inertia. 

The rigid body Lagrangian is the kinetic energy of the system 

L{9.9) = \{9-'9.J{g-^)) = \{^.J0) = m, (3.9) 
where ^ = 9^^9 G so(3) and (•, •) is the pairing between the Lie group and its dual. 



Poisson Structures and Casimir Functions. The Lie algebra dual 5o(3)* has 
a well-known Lie-Poisson structure with a Casimir C^(,(^y{^) = Tr(/i^), where 
/i G so(3)*. Upon identification with 'S? , its generic symplectic leaves become 
concentric spheres with Kirilov-Kostant symplectic form being proportional to the 
area form. If y denotes coordinates on = so(3)*, then the above Casimir function 
is given by CsoO.Y{y) = W vf- 

Following Section 5 of |MPeS 9£ ] on discrete Euler-Poincare reduction, we obtain 
the reduced form of the Moser-Veselov Lagrangian on the group SO (3) given by 

£(/) - Tr(/A), 

where / G SO (3) and SO (3) is embedded into the linear space 0t(3). Then, the 
Legendre transform Ft takes the form 

Fi{f) = L) o MU) = skew(/A) = /A - : S0(3) ^ so(3)*, 

where the constraint that / be in S0(3) has been enforced. The pull-back of C^o(i)* 
under Ff defines a Casimir function on the group, which up to a constant term 
and a sign, is given by 

Cso(3)(/) - Tr(/A/A) / G S0(3). (3.10) 

Its symplectic leaves constitute the invariant manifolds of the reduced discrete dy- 
namics corresponding to the Lagrangian (3.9). 

Analogously, one can define a Poisson structure on the Lie algebra so(3) using 
the duality between Lie-Poisson and Euler-Poincare reduced systems on so(3)* and 
5o(3), respectively. The Lagrangian (|3.9| ) defines the Legendre transformations Fl 

from so(3) to so(3)* given by /i = ^ = J(i). Then, the pull-back by Fl* defines 

a Casimir function on so(3): 

c.o(3)(e) = Fl* ■ c,„(3).(c) = ((j(o,^(e))), 

where the metric on the dual is induced by the one on the algebra, i.e. by the 
symmetric positive definite operator J. If x denotes coordinates on R'^ = so(3), 
then the above Casimir function is given by Cso(3)(a;) = 11^7(0;) |p. Thus, the corre- 
sponding symplectic leaves are ellipsoids of . They do not coincide with adjoint 
orbits, which are spheres in R'^. The dynamic orbits are obtained by intersecting 
these ellipsoids, determined by , with the energy ellipsoids, determined by J. 
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